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Constrained Optimum Path (COP) problems appear in many real-life applications, especially on 
communication networks. Some of these problems have been considered and solved by specific 
techniques which are usually difficult to extend. In this paper, we introduce a novel local search 
modeling for solving some COPs by local search. The modeling features the compositionality, mod- 
ularity, reuse and strengthens the benefits of Constrained-Based Local Search (TJ. We also apply the 
modeling to the edge-disjoint paths problem (EDP). Computational results show the significance of 
the approach. 

1 Introduction 

Constrained Optimum Path (COP) problems, where optimum path from origin to destination satisfying 
additional constraints must be found, appear in many real-life applications, especially on communica- 
tion and transportation networks. These problems have attracted considerable attention from different 
research communities: operations research, telecommunications because of its various applications (see 
El and the references therein). 

Most of COP problems are NP-hard. Some instances can be solved efficiently by specific techniques, 
for instance, branch and bound using a Lagrangian-based bound [3], vertex-labeling H, etc. These 
techniques seem to be sophisticated and depend on particular constraints and objective functions to be 
optimized. Moreover, they are difficult to extend, for example, when we face with generalized problems 
where more constraints are required to be satisfied. We propose in this paper a novel approach for 
modeling and solving some COP problems by local search where the desired paths are elementary (i.e. 
no repeated nodes). The objective of this work is to extend the LS (Graph & Tree) framework Q 
by the design and the implementation of abstractions allowing to easily model and solve some COP 
problems using this approach. The computational model features compositionality, modularity, reuse, 
and strengthens the local search modeling benefits of Constraint-Based Local Search [ 1 ] which provides 
the separation of concerns. 

The proposed computational model has been applied to the EDP problem. Computational results 
show the significance of the approach. 

2 Modeling paths with spanning trees 

We introduce in this section the proposed approach for modeling COP problems with spanning trees. 
We first give somes definitions and notations over graphs. Our framework considers both directed and 
undirected graphs but for simplicity, we consider in this presentation only undirected graphs. 

Given an undirected graph g, we denote V(g), E(g) respectively the set of nodes and the set of edges 
of g. 

Y. Deville and C. Solnon (Eds): Sixth International Workshop 

on Local Search Techniques in Constraint Satisfaction (LSCS'09). © Q. D. Pham, Y. Deville, and P. van Hentenryck 




6 



A Local Search Modeling for Constrained Optimum Paths Problems 



In this paper, we only consider elementary paths, henceforth we use the word "path" instead of 
"elementary path" if there is no ambiguity. A graph is connected if and only if there exists a path from u 
to v, Vm, v G V(g). A tree is an undirected connected graph containing no cycles. A spanning tree tr of an 
undirected connected graph g is a tree spanning all the nodes of g: V(tr) = V(g) and E(tr) C E(g). A 
tree tr is called a rooted tree at r if the node r has been designated the root. Each edge of tr is implicitly 
oriented towards the root. If the edge (w, v) is oriented from u to v, we call v the father of u on tr. 

The key decision design is inspired from the following observation: Given a rooted tree tr whose 
root is t, the path from a given node s to t on tr is unique. An update of tr will generate a new rooted tree 
which may induce a new path from s to t on this tree. 

Given an undirected graph g and a node r G V(g), VarRootedSpanningTree(g,r) (also called rooted 
spanning tree variable) is a concept representing a dynamic spanning tree of g. The spanning tree is 
rooted at r. 

In order to model a COP problem in an undirected graph g in which the source and the target of 
the desired path are respectively s,t G V(g), we use VarRootedSpanningTree(g,s,t) which is VarRoot- 
edSpanningTree(g ,t) with a node s designated as the source node. Each instance tr of VarRootedSpan- 
ningTree(g,s,t) specifies a unique path from s to t on g. Henceforth we use s to denote the source node 
of any rooted spanning tree of the given graph and the path from s to the root of a rooted spanning tree 
tr on tr is called the path induced by tr if there is no ambiguity. 

The main avantage of using rooted spanning tree for modeling paths instead of using explicit paths 
representation (i.e. a sequence of nodes) is the simplification of neighborhood computation. The tree 
structure constains rich information that induces directly path structure from a node s to the root. A 
simple update over that tree (i.e. an edge replacement which is detailed in Section [3} will induce a new 
path from s to the root. 

3 Neighborhood 

Given an instance tr of VarRootedSpanningTree(g,s,t), we show how to change tr in order to generate a 
new rooted spanning tree tr 1 of g which induces a new path from s to t on g. 

Given an undirected graph g, an instance tr of VarRootedSpanningTree(g,s,t), an edge e = (w,v) 
such that e € E(g) \ E(tr) is called replacing edge of tr. We denote rpl(tr) the set of replacing edges 
of tr. Given e G rpl(tr), an edge e' that belongs to the path between two endpoints of e on tr is called 
replacable edge of e. We denote rpl(tr,e) the set of replacable edges of e. Intuitionally, a replacing edge 
e is an edge that is not in the tree tr but that can be added to tr (this edge insertion creates a cycle C when 
we ignore orientations of edges of tr), and all edges of this cycle except e are replacable edges of e. 

Given an undirected graph g, an instance tr of VarRootedSpanningTree(g,s,t), e and e' are respec- 
tively replacing edge of tr and replacable edge of e, we define the following edge replacement action: 

1. Insert the edge e = (u,v) to tr. This creates an undirected graph g with a cycle C containing the 
edge e'. 

2. Remove e' from g. 

After taking above edge replacement action, we obtain a new rooted spaning tree tr' of g. We denote 
tr' = rep(tr,e' ,e). The neighborhood of tr is 

N(tr) = {tr = rep(tr,e' ,e) \ e G rpl(tr),e' G rpl(tr,e)} 

It is easy to observe that two different spanning trees tr\ and tr2 rooted at t of an undirected graph g 
may induce the same path from s to t on g (s,t G V(g)). The neighborhood N(tr) must then be reduced 
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a. current tree tr b. tr' = rep(tr, (7, 11), (8, 10)) 



Figure 1 : Example of basic move 





© © 

a. current tree tr b. tr' = rep(tr, (7, 11), (8, 10), (3,4), (1,5)) 



Figure 2: Example of complex move 



such that the new tree induces a new path from s to t. This reduction is described in The action 
rep(tr,eo,ei) is called a basic move. Figure 1 gives an example of basic move. 

It is possible to consider more complex moves by applying a set of independent basic moves. Two 
basic moves are independent if the execution of the first one does not affect the second one and vice versa. 
The execution order of these basic moves does not affect the final result. Figure 2 gives an example of 
complex move. 



4 COMET implementation 

We extend the LS (Graph & Tree) framework by implementing some Graphlnvariants, GraphCon- 
straints and GraphObjectives (see for more detail) for modeling and solving some COP problems. 
Graphlnvariant is a concept representing objects which maintain some properties of a dynamic graphT 



dynamic graph is a graph that can be changed e.g., by the removal or the insertion of vertices, edges. 
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1. LSGraphSolver ls(); 

2. VarRootedSpanningTree tr(ls,g,s,t) ; 

3. Pref erredReplacingEdges pref Replacing(tr) ; 

4. Pref erredReplacableEdges pref Replacable (tr) ; 

5. forall(ei in pref Replacing. getSet () ) 

6. forall(eo in pref Replacable. getSet(ei)) 

7. d = C.getReplaceEdgeDelta(tr,eo,ei) ; 

Figure 3: Exploring the basic neighborhood 



(for instance, the sum of weights of all the edges of a graph, the diameter of a tree, etc.). GraphConstraint 
and GraphObjective are concepts describing differentiable objects which maintain some properties (for 
instance, the number of violations of a constraint or the value of an objective function) of a dynamic 
graph. The main feature of GraphConstraint and GraphObjective is the interface, allowing to query 
the impact of local moves (modification of the dynamic graph) on these properties. Some GraphCon- 
straints and GraphObjectives have been designed and implemented over VarRootedSpanningTree. For 
instance, PathCostOnEdges(tr ,k) |^]is an abstraction representing the total weight accumulated along 
the path induced by tr, MinEdgeCost (tr ,k) , MaxEdgeCost(tr,k) are abstractions representing the 
minimal and maximal weight of edges on the path induced by tr. NodesVisited(tr , S) is an abstrac- 
tion representing the number of nodes of S visited by the path induced by tr. These abstractions are 
example of GraphObjectives which are fundamental when modeling Constrained Optimum Path prob- 
lems. For example, in QoS, we consider shortest path from an origin to a destination with constraints over 
bandwidth which is defined to be the minimum weight of edges on the specified path. The framework 
provides flexibility for modeling various Constrained Optimum Path problems. We can easily combine 
(with +,-,* operator) and state basic constraints (for instance, <=,>=,==) over these abstractions. 
PathEdgeDis joint (tr) is another GraphConstraint which is defined over an array of paths (the i 
path is induced by tr [i] ) and specifies that these paths are mutually edge-disjoint. 

In order to illustrate the modeling, we give a snippet (see Figure 3) which explores the basic neigh- 
borhood. Line 1 initializes a LSGraphSolver object Is which manages all the VarGraph, VarTree, 
Graphlnviants, GraphConstraint?, and GraphObjectives and relations (dependency graph) between these 
objects. Line 2 declares and initializes randomly a VarRootedSpanningTree tr rooted at t of the input 
graph g which represents the path from the source node s to t. pref Replacing and pref Replacable 
are Graphlnvariants which maintain the set of preferred replacing edge^jand preferred replacable edges 
(lines 3-4). Lines 5-7 explore the basic neighborhood and evaluate the quality of moves with respect to a 
GraphConstraint C. The getReplaceEdgeDelta (line 7) method returns the variation of the number of 
violations of C when the preferred replacable edge eo is replaced by the preferred replacing edge ei on 
tr. 



5 Application: The EDP problem 

Given an undirected graph G = (V,E), and a set T = {< Sj,tj >| st ^ G V} representing a list of 
commodities (jjr = k). EDP consists of finding a maximal cardinality set of mutually edge-disjoint paths 

2 tr is a VarRootedSpanningTree, k is the index of the considered weight on edges. 

3 edges used for the edge replacements which return new trees inducing new paths from s to t. 
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void stateModel{ 

1. LSGraphSolver ls(); 

2. VarRootedSpanningTree tr[i in 1 . .k] (Is ,g, s [i] ,t [i] ) ; 

3. PathEdgeDis joint ed(tr) ; 

4 . Is . close () ; 

5. } 

Figure 4: The Modeling for the EDP problem 

from Sj to t, on G (< Si,t\ >£ T). In Q, a Multi-start Simple Greedy algorithm is presented as well 
as an ACO algorithm. The ACO is known to be state-of-the-art for this problem. We propose a local 
search algorithm using the modeling approach. The model is given in Figure 4 where line 2 initializes 
an array of k VarRootedSpanningTrees representing k paths between commodities. The edge-disjoint 
constraint ed is which is defined over paths from s [i] to t [i] on tr [i] (i = 1 , 2 , . . . , k) stated 
in line 3. In [7 ], the following criterion is introduced which quantifies the degree of non-disjointness of 
a solution 5 = {P\,Pz, ...fit} (Pj is a P atri from Sj to tj): 

C(S) = £(max{0, £ p'(S,e) - 1}) 

eeE PjeS 

where pi(S,e) = 1, if e G Pj 6 5 and pj(S,e) = 0, otherwise. 

The number of violations of the PathEdgeDis joint (Pi, P2, ...,P/t) constraint in the framework is de- 
fined to be C({Pi,Pz, ...,-Pfc}) and the proposed local search algorithm tries to minimize this criterion. 
From a solution which is normally a set of k non-disjoint, a feasible solution to the EDP problem can 
be extracted by iteratively removing the path which has most edges in common with other paths until 
all remaining paths are mutually edge-disjoint as suggested in [7]. In our local search model, we extend 
that idea by taking a simple greedy algorithm over the remaining paths after that extraction procedure in 
hope of improving the number of edge-disjoint paths. 

The main idea for the search is to try different moves to get a first improvement: 1-move or 2-move 
over one VarRootedSpanningTree, two 1-moves at hand over two VarRootedSpanningTrees. 

For the experimentation, we re-implemented in COMET the Multi-start Greedy Algorithm (MSGA) 
and the ACO (the extended version) algorithm which are described in Q, and compare them with our 
local search model. The instances experimented (graphs including commodities) in the original paper 
171 is not available at the moment (except some graphs). We base on the instances generation description 
in the paper (7J and generate new instances as follows. We take 4 graphs from Q. For each graph, 
we generate randomly different sets of commodities with different sizes depending on the size of the 
graph: for each graph of size n, we generate randomly 20 instances with 0.10*n, 0.25*?i and 0A0*n 
commodities. In total, we have 240 problems instances. Due to the huge complexity of the problem, we 
execute each problem instance once with a time limit of 30 minutes for each execution. Experimental 
results are shown in Table 1. The time window for the MSGA and the ACO algorithms are also 30 
minutes. The Table reports the average values of the objective function of the best solutions found and 
the average of times for obtaining these best solutions of 20 instances (a graph G = (V,E) and a set of 
r*\V \ commodities, r = 0.10,0.25,0.40). The Table shows that our local search model gives competitive 
results in comparison with the MSGA and the ACO algorithms. In comparison with the MSGA, our local 
search model find better solutions in 217/240 instances while the MSGA find better solutions in 4/240 
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Table 1 : Experimental results of EDP problem 



instances. On the other hand, in comparison with the ACO model, our local search model find better 
solutions in 144/240 instances while the ACO model find better solutions in 11/240 instances. 



6 Conclusion and future work 

We introduce in this paper a novel local search modeling for Constrained Optimum Path problems on 
graphs. The objective here is to give a high-level modeling framework for implementing some COP 
problems which strengthens the benefits modeling of CBLS and features compositionality, modularity 
and reuse. The modeling provides a clean seperation of concerns: The modeling component and the 
search component are independent. On one hand, it is easy to add new constraints and to modify or re- 
move existing ones, without having to worry about the global effect of these changes. On the other hand, 
programmers can experiment with different heuristic and metaheuristics without affecting the problem 
modeling. The modeling is based on VarRootedSpanningTree(g,s,t) concept inspiring the observation 
that each tree induces a unique path between two specified nodes and a update over this tree satisfying 
additional contraints generates a new tree which induces a new path between these nodes. Various neigh- 
borhoods for have been defined. The abstraction is implemented by extending the LS (Graph & Tree) 
in COMET. The modeling has been experimented on the resources constrained shortest path problem and 
the edge-disjoint paths problem which show the significance of the framework. 
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